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Context. The chemical composition of the gas in galaxies versus cosmic time provides a very important tool for understanding galaxy 
evolution. Although there are many studies at high redshift, they are rather scarce at lower redshifts. However, low redshift studies 
can provide important clues about the evolution of galaxies, furnishing the required link between local and high redshift universe. In 
this work we focus on the metallicity of the gas of star-forming galaxies at low redshift, looking for signs of chemical evolution. 
Aims. To analyze the metallicity contents star-forming galaxies of similar luminosities and masses at different redshifts. With this 
purpose, we present a study of the metallicity of relatively massive (log(M„ ar /M ) > 10.5) star forming galaxies from SDSS-DR5 
(Sloan Digital Sky Survey-Data Release 5), using different redshift intervals from 0.04 to 0.4. 

Methods. We used data processed with the STARLIGHT spectral synthesis code, correcting the fluxes for dust extinction, estimating 
metallicities using the R 2 i method, and segregating the samples with respect to the value of the [Nil] /16583/[0 n] A3727 line ratio in 
order to break the R 2 i degeneracy selecting the upper branch. We analyze the luminosity and mass-metallicity relations, and the effect 
of the Sloan fiber diameter looking for possible biases. 

Results. By dividing our redshift samples in intervals of similar magnitude and comparing them, significant signs of metallicity 
evolution are found. Metallicity correlates inversely with redshift: from redshift to 0.4 a decrement of ~0.1 dex in 12+log(0/H) is 
. found. 

■ Key words, galaxies: abundances - galaxies: evolution - galaxies: starburst 
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r*"** ' 1 . Introduction method to estimate metallicities in galaxies is known as the "T e 

O I method" (Pagel et al. 1992; Skillman & Kennicutt 1993), which 

Ch , Determination of the chemical composition of the gas and stars consists on measuring the ratio of the [Q m] ^4353 auroral line t0 

O ■ m galaxies versus cosmic time provides a very important tool a lower excitation line such as [Q m] ^ 5007 Assuming a classi- 
^ ; for understanding galaxy evolution, due to its important impact cal H n region model> Ms ratio provides an estimate of the elec . 
. - . on fields such as stellar evolution and nucleosynthesis, gas en- tfon temperature of the gas , which is then conver ted into metal- 
^ ■ nchment processes, and the primary or secondary nature of the mty (0sterbrock 1987) . However, [Om] ,14363 is too weak to 
H i different chemical species. Historically, the main observational be easily observedi not only in metal ^ch, but even in metal 
«? , evidence suggestive of chemical evolution of galaxies, is pro- galaxies, Z < 0.5 Z (log(0/H) + 12 < 8.6), and accord- 
vided by the observation of different chemical compositions of ing tQ Kobulnicky et al . (1999X for low-metallicity galaxies, 
stars of different ages of the Milky Way and its environment (see the [Q m] M363 diagnostic system atically underestimates the 
the reviews by Audouze & Tinsley 1976, Wheelet et al. 1989, global oxygen abundance . Als0; the sam e authors show that for 
Wilson & Matteucci 1992, Mc William 1997). massive, metal-rich galaxies, empirical calibrations using strong 
The study of the evolution of the metal enrichment in galax- emission-line ratios can be reliable indicators of the overall oxy- 
ies is based mainly in two methods. One is based on the detection gen abundance in H n regions 
of absorption lines in QSO spectra produced by the neutral in- 
terstellar medium (ISM) of galaxies in the line-of-sight of the For these reasons, theoretical metallicity calibrations of 
QSO (Prochaska et al. 2003), while the other uses emission lines strong-line ratios using photoionization models are used instead 
of the warm ISM (H 11 regions) detected in the integrated galaxy for determining abundances of high metallicity star-forming 
spectra. galaxies, such as: [Nn] /i6584/[Oii] A3121 (Kewley & Dopita 
Optical emission lines in galaxies have been widely used 2002) and the R23 ratio, introduced by Pagel et al. (1979). The 
to estimate abundances in extragalactic Hn regions (e.g. Aller first one provides an excellent abundance diagnostic for Z > 
1942; Searle 1971; Pagel 1986; Shields 1990, among others). 0.5 Z (log(0/H) + 12 > 8.6), because N + and + have simi- 
Among the different methods developed to estimate metallici- lar ionization potentials, and this ratio is almost independent of 
ties, we can distinguish between theoretical models, empirical ionization parameter. However it cannot be used at lower abun- 
calibrations, or a combination of both (for a review see, e.g., dances (Z < 0.5 Z ), where the metallicity dependence of the 
Kewley & Dopita 2002; Kewley & Ellison 2008). The direct [Nn] ^6584/[On] ,0727 ratio is lost because nitrogen (like oxy- 
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gen) is predominantly a primary nucleosynthesis element in this 
metallicity range (Kewley & Dopita 2002). The /?23 method, is a 
widely used and well calibrated method (see for example Alloin 
et al. 1979; Edmunds & Pagel 1984; McCall et al. 1985; Dopita 
& Evans 1986; McGaugh 1991; Zaritsky et al. 1994; Kewley & 
Dopita 2002; Kobulnicky & Kewley 2004; Tremonti et al. 2004, 
hereafter T04; Liang et al. 2006, hereafter L06). However, it 
has the disadvantages of being double-valued as a function of 
12+log(0/H), and that depends on the ionization parameter, par- 
ticularly for Z < 0.5 Z , being less sensitive to metallicity in this 
range. 

Alternatively, when the direct method can not be used, em- 
pirical calibrations can be obtained by fitting the relationship 
between direct T e metallicities and strong-line ratios as well. 
Typical empirical calibrations are: the /?23 ratio (Pilyugin 2001; 
Pilyugin & Thuan 2005; Liang et al. 2007), from which Pilyugin 
(2001) derived an empirical calibration based on T e metallicities 
for a sample of Hn regions, the [Nn] /16583/ Ha ratio (Pettini 
& Pagel 2004, hereafter PP04), and the ([Om] /l5007/H J S)/([Nii] 
^6583/Ha) ratio (03N2), (PP04). Although the latest method is 
of little use when 03N2 > 2, at lower values the relation is rela- 
tively tight and linear (PP04). 

Finally, as an example of a combined calibration, we have 
the N2 = [Nil] ,16583/ Ha method (Denicolo et al. 2002), which 
follows a linear relation with log(0/H) that holds for approx- 
imately from l/50th to twice the Solar value. This method is 
based on a fit to the relationship between the T e metallicities 
and the [Nn] /16583/Ha ratio, from which some have metallici- 
ties derived using the T e method, and the remaining metallicities 
were estimated using either the theoretical /?23 or an empirical 
method. 

Nevertheless, comparisons among the metallicities estimated 
using different theoretical and empirical methods reveal large 
discrepancies (e.g., Pilyugin 2001, Bresolin et al. 2004, Garnett 
et al. 2004), with theoretical calibrations favouring higher metal- 
licity values than those obtained using electron temperature es- 
timations. 

In the field of metallicity evolution versus cosmic time, there 
exist many studies, both theoretical and observational. Among 
the models we have, for example, that of Buat et al. (2008) and 
Kobayashi et al. (2007), who derived models of metallicity as a 
function of z, which show a progressive increase in metallicity 
with time, even at low redshift. Savaglio et al. (2005), devel- 
oped an empirical model of metallicity evolution based on ob- 
servations, in which the metallicity at z < 1 is an interpolation of 
that at higher redshifts. Also, Brooks et al. (2007) and Finlator 
& Dave (2008), among others, derived cosmological models of 
the mass-metallicity relation. The cosmic metal enrichment is 
attributed to a higher past volume-averaged star formation rate 
(see for example Madau et al. 1996; Lilli et al. 1996; Flores et 
al. 1999). 

In addition, the metallicity and masses of galaxies are 
strongly correlated, with massive galaxies showing higher metal- 
licities than less massive galaxies. This mass-metallicity (M-Z) 
relation has been intensively studied (Skillman et al. 1989; 
Brodie & Huchra 1991; Zaritsky et al. 1994; Richer & McCall 
1995; Garnett et al. 1997; Pilyugin & Ferrini 2000, among oth- 
ers), and it is well established in the local universe (z ~ 0.1) 
by the work of T04 using SDSS data. Regarding to this evolu- 
tion of the mass-metallicity relation of star-forming galaxies at 
high redshift, Erb et al. (2006) found that star-forming galax- 
ies at redshift ~ 2 have 0.3 dex fainter metallicities. Similarly, 
Maiolino et al. (2008) found evolution at z ~ 3.5, which appear 
to be much stronger than the one observed at lower redshifts, 



suggesting that this redshift corresponds to an epoch of major 
activity in terms of star formation and metal enrichment. At in- 
termediate redshifts (1 < z < 2), there are several important stud- 
ies of the evolution of the chemical composition of the gas, such 
as the ones by Maier et al. (2006), Perez-Montero et al. (2009), 
and Liu et al. (2008), the last one found that the zero point of the 
M-Z relation evolves with redshift, in the sense that galaxies at 
fixed stellar mass become more metal-rich at lower redshift. 

Among the studies at z < 1, usually based on small samples, 
Savaglio et al. (2005) have investigated the mass-metallicity re- 
lations using galaxies at 0.4 < z < 1, finding that metallicity is 
lower at higher redshift, for the same stellar mass, by ~ 0.15 
dex. Also, Maier at al. (2005), from a sample of 30 galaxies with 
0.47 < z < 0.92, found that one-third have metallicities lower 
than those of local galaxies with similar luminosities and star 
formation rates. Consistently, Hammer et al. (2005) and Liang 
et al. (2006) found that at z ~ 0.7, emission line galaxies were 
poorer in metals than present-day spirals by 0.3 dex. However, 
Kobulnicky & Kewley (2004) report a smaller variation of 0.14 
dex for < z < 1. This difference could be attributed to the 
fact that the last authors used equivalent widths and standard 
underlying stellar absorption, rather than high quality calibrated 
spectra for measuring the Balmer absorption. Lilly et al. (2003), 
from a sample of 66 star forming galaxies with 0.47 < z < 0.92, 
found a smaller variation in metallicity of ~ 0.08 dex compared 
with the metallicity observed locally, showing only modest evo- 
lutionary effects. On the contrary, Carollo & Lilly (2001), from 
emission-line ratios of 15 galaxies in a range of 0.5 < z < 1, 
found that their metallicities appear to be remarkably similar to 
those of local galaxies selected with the same criteria. A simi- 
lar result, consistent with no significant evolution, was found for 
the luminosity-metallicity relation by Lamareille et al. (2006), 
comparing 131 intermediate redshift star-forming galaxies (0.2 
< z < 1, split in 0.2 redshift bins). However, a recent study of 
Lamareille et al. (2009), focussed on the evolution of the M-Z 
relation up to z ~ 0.9, suggesting that the M-Z relation is 
flatter at higher redshifts. At z ~ 0.77, galaxies of 10 9 4 solar 
masses have -0.18 dex lower metallicities than galaxies of sim- 
ilar masses in the local universe, while galaxies of 10 10 2 solar 
masses have -0.28 dex lower metallicities. 

These discrepancies point to a need for studying the lower 
redshift galaxy samples, to ascertain whether or not at such low 
redshifts (i.e. an age of 8.4 Gyr for z ~ 0.5, using a concordance 
A-CDM cosmology, H =70, O ffl = 0.3 and Q. A = 0.7; Spergel 
et al. 2003) there exist evidences for metallicity evolution, as 
also to serve as calibrator of higher redshift studies. However, 
to be able to compare different redshift samples, it is advisable 
to use the same method for estimating metallicities, since, as 
explained above, theoretical and empirical calibrations generate 
discrepancies in the metallicity estimates depending on the used 
method. 

The SDSS database provides an excellent opportunity for 
extending these studies down to z ~0.4, in order to explore a 
possible evolution of metallicity at low-redshift, but using large 
samples, thus deriving more statistically significant results. In 
this paper, we extend the study presented in our previous article 
(Lara-Lopez et al. 2009, hereafter L09), from 207 to more than 
12000 galaxies, spanning more luminosity intervals in redshift 
bins of 0. 1 from ~0 to 0.4, and adding analyses of the mass and 
luminosity-metallicity relations, as well as about the origin of 
nitrogen in our galaxies. 

This paper is structured as follows: in Sect. 2 we give a de- 
tailed description of the used data, in Sect. 3 we describe the 
metallicity estimates, the [Nn] /16583/[Oii] A3121 diagram and 
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its metallicity distribution, in Sect 4 we investigate the origin 
of nitrogen of our galaxies, whereas in Sect. 5 we discuss our 
results taking into account the possible biases of our samples. 
Finally, conclusions are given in Sect. 6. 

2. Sample selection 

We analyzed the properties of a selected sample of emis- 
sion lines galaxies from SDSS-DR5 (Adelman-McCarthy et al. 
2007). Targets were observed using a 2.5 m telescope located at 
Apache Point Observatory (Gunn et al. 2006). The SDSS spec- 
tra were obtained through 3 arcsec diameter fibres, covering a 
wavelength range of 3800-9200 A, and with a mean spectral 
resolution A/AA ~ 1800. The SDSS-DR5 spectroscopy database 
contains spectra for ~ 10 6 objects over ~ 5700 deg 2 . Further 
technical details can be found in Stoughton et al. (2002). 

We used the SDSS-DR5 spectra from STARLIGHT 
Database^], which were processed through the STARLIGHT 
spectral synthesis code, developed by Cid Fernandes and col- 
leagues (Cid Fernandes et al. 2005, 2007; Mateus et al. 2006; 
Asari et al. 2007). From them, we obtained the emission lines 
fluxes measurements of our samples from the contiuum sub- 
tracted spectra. For each emission line, STARLIGHT code re- 
turns the rest frame flux and its associated equivalent width, 
linewidth, velocity displacement relative to the rest-frame wave- 
length and the S/N of the fit. In the case of Balmer lines, 
STARLIGHT code corrects for underlying stellar absorption us- 
ing synthetic spectra obtained by fitting an observed spectrum 
with a combination of 150 simple stellar populations (SSPs) 
from the evolutionary synthesis models of Bruzual & Chariot 
(2003), computed using a Chabrier (2003) initial mass function, 
"Padova 1994" evolutionary tracks (Alongi et al. 1993; Bressan 
et al. 1993; Fagotto et al. 1994a, b; Girardi et al. 1996), and 
STELIB library (Le Borgne et al. 2003). The 150 base elements 
span 25 ages between 1 Myr and 18 Gyr, and six metallicities, 
from Z =0.005 to 2.5 Z ; for more details see Mateus et al. 
(2006). 

Our objective is to study the chemical evolution of emission- 
line galaxies by comparing galaxies at different redshifts in small 
and equal ranges of luminosities. To this aim, our sample is di- 
vided in redshift intervals of 0.1 from z ~ to 0.4. The pur- 
pose for selecting small ranges in luminosity is to alleviate the 
problem of the magnitude completeness, since the Sloan sam- 
ple is only complete in the magnitude range 14.5 < m r < 17.7 
(e.g., Asari et al. 2007), and then becomes incomplete at redshift 
above z > 0. 1 (e.g., Kewley et al. 2006). However, this procedure 
limits the study to the more luminous galaxies, since they are the 
ones detected at any redshift interval. 

With this aim, our initial sample was divided in the following 
redshift intervals: z =(0.04-0.1), z\ =(0.1-0.2), z 2 =(0.2-0.3) 
and Z3=(0.3-0.4). To ensure covering > 20% of the light, we 
selected galaxies for the zo sub-sample with z > 0.04 , as re- 
comended by Kewley et al. (2005). This preliminary selection 
give us 197967 galaxies for zo, 226012 for zi, 42205 for zi, and 
38305 for the Z3 sub-samples. Absolute magnitudes were both 
K and Galactic extinction corrected, by using the code provided 
by Blanton et al. (2003), and the maps of Schlegel et al. (1998), 
respectively, as provided by the STARLIGHT team. 

In order to determine the luminosity-intervals at each red- 
shift, we proceed to estimate the completeness of the absolute 
Petrosian r magnitude for each redshift-interval. We obtained 
our luminosity-complete samples selecting all galaxies brighter 

1 http://www.starlight.ufsc.br 



than the median for the larger redshift interval to be com- 
pared with, as shown in Fig. 1. According to that, we proceed 
constructing three luminosity-samples (hereafter 1-samples): 
the 1-sample a was constructed by taking the interval of the 
luminosity-completeness of z\ (-23.8 < M r < -21.7), this al- 
lows comparing redshift-intervals zo, and z\, in 1-sample b we 
take the luminosity-completeness of zi (-24.8 < M r < -22.9), 
for comparing redshift-intervals zi, and z\ ; and in 1-sample c we 
take the luminosity-completeness of zj, (-24.8 < M r < -23.1), 
to compare redshift-intervals zs, Zi, and z\, as shown in Table 1 
and Fig. 1. As resulting from this figure, it was not possible to 
introduce the zo redshift-sample for comparison with 1-samples 
b and c, due to the small number of galaxies at their luminosity 
completeness. 



L-sample 


redshift 




M r 


redshift-samples 




completeness 






compared 


a 


Zl 


(" 


•23.8-21.7) 


Zl, zo 


b 


Z2 


(" 


■24.8,-22.9) 


Zl, Zl 


c 


Z3 


(" 


•24.8,-23.1) 


Zi, Z2, Zl 



Table 1. Luminosities-samples studied with their respective lu- 
minosity interval in absolute Petrosian r-magnitude, and the 
redshift-samples that can be compared. 



The resulting samples contain, for the 1-sample b: 2352 
galaxies for zi and 1386 for z\\ for 1-sample a: 31111 for zi and 
9288 forzo- The 1-sample c, was already studied in in L09. 

From these samples we only consider galaxies whose spectra 
show in emission the Ha, H/3, [Nn] /16583, [On] A3121, [Ora] 
,14959 and [Om] A5001 lines, with Ha, HyS, [Nn] ,16583 and 
[O n] A3121 signal-to-noise ratios higher than 3cr. 

Finally, we selected star-forming galaxies following the 
criteria given by Kauffmann et al. (2003) in the Baldwin 
et al. (1981) diagram: log[Ora] A5001/H/3 < 0.61/{log([Nn] 
/i6583/Ha)-0.05} + 1.3, used for example by Veilleux & 
Osterbrock (1987), Kewley et al. (2001, 2006), Stasihska et al. 
(2006), among others. After all these selections, the number of 
galaxies for each sample is, for the 1-sample b: 335 galaxies for 
Zz, and 148 for zi, and for the 1-sample a: 10477 for zi, and 1577 
for zo- 



2. 1 . Dust extinction 

Since Balmer lines are already corrected for underlying stellar 
absorption by the STARLIGHT code, it only remains correcting 
for dust extinction. Our extinction correction was derived using 
the Balmer decrements in order to obtain the reddening coeffi- 
cient C(HjS). Assuming case-B recombination with a density of 
100 cirT 3 and a temperature of 10 4 K, the predicted ratio (unaf- 
fected by reddening or absorption) of Ha/Hfi is 2.86 (Osterbrock 
1989), and the coefficient is given by: 



where F(A) and 1(A) are the observed and the theoretical fluxes, 
respectively, and f(A) is the reddening curve normalized to H/3 
using the Cardelli et al (1989) law, with R v = A,,/E(B - V)=3.1 
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3. Metallicity estimates and evolution 




Fig. 1. Absolute Petrosian r-magnitude histograms. From top to 
bottom, histograms for the 23, zi, Z\, and zo sub-samples, re- 
spectively. Dot-dashed line represents the completeness for the 
Z3 sample, solid line shows the completeness for the zi sample, 
and dashed lines the completeness for the z\ sample. All of them 
extended along the 1-samples a, b and c, where N is the number 
of galaxies, respectively. 

Once we have obtained the reddening coefficient for each 
galaxy of our samples, we proceed to estimate the corrected 
fluxes using F corr (A)= F obs (A)l0 OAA \ with, 

A w = 2.5 C(H(3) 
A [N n]i6583 = 1.747 C(H/3) 
A Ha = 1.758 C(Hp) 
A [O iii]^5007 = 2.403 C(H(3) 
A[omn4959 = 2 - 43 3 C(flJS) 
A [0 ii],i3727 = 3.303 C(Hj3), 

as calculated from a prescription given by Cardelli et al. (1989) 



We estimate metallicities using the R23 relation, introduced by 
Pageletal. (1979), 

fl23=([On] ^3727+[Om] ,1,14959, 5007)/H/3, 

adopting the calibration given by Tremonti et al. (2004), 

12 + log(0/H) = 9.185 - 0.313x - 0.264x 2 - 0.321x 3 , (1) 

where x = log R23 ■ 

However, this calibration is valid only for the upper branch 
of the double-valued R23 abundance relation, and additional line 
ratios, such as [Nn] /i6583/[On]/13727, are required to break this 
degeneracy. Since the upper and lower branches of the R23 cali- 
bration bifurcate at log([Nn]/[On]) ~ -1.2 for the SDSS galax- 
ies (Kewley & Ellison 2008), which corresponds to a metal- 
licity of 12+log(0/H) ~ 8.4, we further select galaxies having 
12+log(0/H) > 8.4 and log([Nii]/[On]) > -1.2, corresponding 
to the upper R23 branch. Applying this final discrimination, we 
end for 1-sample b with 331 galaxies for 12 and 146 for z\, and 
for 1-sample a with 10434 galaxies for z\ and 1576 for zo- These 
are the samples that will be analyzed in this paper. 

From these data, we derived the abundance-sensitive diag- 
nostic diagram [Nn] ,16583/[On] ,13727 vs. 12+log(0/H), rep- 
resented in Fig. 2. This diagram has also been used, for ex- 
ample, by Kewley & Dopita (2002), Nagao et al. (2006), and 
L06, among other metallicity-sensitive emission-line ratios, like 
log([Nn] ,16583/Ha), log([Ora] i5007/H/3)/[N 11] i6583/Ha-), 
and log([Om] ,14959, 5007/H/3). We selected this specific dia- 
gram due to its low scatter and to the physical additional in- 
formation that provides. The advantages of using [Nn] ,16583 
and [O n] ,13727 lines, is that they are not affected by underly- 
ing stellar population absorption, and because this ratio is almost 
independent of the ionization parameter, since N + and + have 
similar ionization potentials. From all the diagnostic diagrams 
cited before, this one presents the lowest scatter, and since both 
axis are sensible to metallicity, possible signs of evolution could 
be easily identified in it. 

For every 1-sample of Fig. 2, we fit both a linear and an order 
three polynomial, obtaining the coefficients shown in Table 2. 





Linear Fits 




Polynomial Fits 


1-sample 


b b\ 


a 


a a i a 2 


a 


9.139 0.575 


0.097 


9.117 0.329 -0.447 


b 


9.130 0.564 


0.091 


9.112 0.408 -0.226 


c 


9.137 0.599 


0.114 


9.114 0.438 -0.183 



Table 2. Coefficients for the 1-samples. For linear fits we assume 
y — bo + b\x, and for polynomial fits y — ao + a\X + a2X 2 , with y 
= 12+log(0/H) and x = log([Nn] ^6583/[On] ,13727). 



To interpret our results, it is important to note that we are 
working with the integrated spectra, and it is well known that 
the metallicity decreases with the distance from the galaxy cen- 
ter (Garnett et al. 1997). As shown by Kewley et al. (2005), and 
L06 using a sample of Jansen et al. (2000), data points from 
nuclear spectra follow the SDSS galaxies nuclear spectra very 
well, but data points from the integrated spectra show lower 
12+log(0/H). Nuclear metallicities exceed metallicities derived 
from integrated spectra by ~0.13 dex on average. 
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Fig. 2. Calibration relation between 12+log(0/H) and log([Nn] 
/16583/[Oii] /B727). From top to bottom, 1-sample a, b, and c. 
Circles, triangles, squares and asterisks, represent galaxies for 
the redshift-sample zo, Zi, zi, and Z3, respectively. The solid line 
represents the best fit of the data using a linear fit, with the 2<x 
dispersion indicated by the short-dashed lines, while the dot- 
dash line shows the order-3 polynomial fit. 



As shown in L09, for the 1-sample c there is a clear decre- 
ment in the Z3 redshift-sample of ~0.1 dex in 12+log(0/H) with 
respect to the z\ and zi redshift-samples. However, 1-samples a 
and b, show only small decrements, indicating that the redshift 
0.3 represents an important epoch in the evolution of galaxies. 
As argued by Kewley et al. (2008), using a single metallicity 
calibration, the difference in relative metallicities should be the 
same using any other metallicity calibration, although the ab- 



Fig. 3. Normalized metallicities histograms for our three 1- 
samples. Dotted line represents the redshift-sample interval zo, 
dashed line the redshift-sample z\ , solid line the redshift-sample 
zi and dot-dash line the redshift-sample zi- 



solute metallicities might differ from one calibration to another. 
Thus, our main result is the relative decrement in metallicity of 
~0.1 dex of Z3 with respect to the zi and zi redshift-samples of 
1-sample c. 

In order to study the behaviour of the metallicities for the dif- 
ferent 1-samples, we proceed to generate a metallicity histogram 
for our three 1-samples as shown in Fig. 3. In the histograms we 
can observe a shift to lower metallicities as redshift increases, 
which is more evident in 1-sample c, as argued in L09. L-samples 
a and b do not show any important variation in metallicity, as 
can be seen in Table 3. As a measure of the dispersion of the his- 
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tograms, we estimated the interval that encompasses 66% of the 
galaxies around the mode of the distribution. With this criterium, 
we find for 1-sample a, an interval of 0.16 dex in 12+log(0/H) 
for zo, and 0.17 dex for z\\ for 1-sample b, 0.17 dex for zu and 
0.16 for zi, and for 1-sample c, 0.26 dex for zi, 0.18 for zi, and 
0.35 for Z3- As can be observed in Table 3, the dispersion of the 
metallicity histograms increases for the samples with the more 
limited number of galaxies, as can be seen in 1-sample c for zi 
andz3. 

For the redshift-samples we can observe a maximum metal- 
licity for zo, followed by a small decrement in metallicity for 
redshift z\ which remains constant for zi, and then an important 
decrement for zj,. As indicated by Carollo et al. (2001), the red- 
shift interval 0.5 < z < 1 represents a transition between the 
high-redshift universe at z > 1 and that seen today. Then evolu- 
tionary effects should be more evident in the galaxy population 
at these redshifts. Nevertheless, in this work we find a significant 
metallicity evolution at redshift 0.3. 

4. Nitrogen and Oxygen abundances 

The primary and secondary origin of nitrogen is of high impor- 
tance to understand the processes inside stars, and the evolution 
of galaxies. Although this is not the aim of this paper, we are 
in a position to investigate the origin of nitrogen in our galaxy 
samples. 

The nuclear mecanism producing nitrogen in stellar interi- 
ors result from the CN cycle of the CNO reactions, which takes 
place in the stellar hydrogen burning layer, with the net result 
that 14 N is synthesized from 12 C and lfi O (Meynet & Maeder 
2002; Pettini et al. 2008). Nitrogen can be of either primary or 
secondary origin. If the oxygen and carbon are produced in the 
star prior to the CNO cycling, then the amount of nitrogen pro- 
duced is said to be primary. If initial amounts of oxygen and 
carbon are incorporated into a star at its formation, and a con- 
stant mass fraction is processed, then the amount of nitrogen 
produced is proportional to the initial heavy-elements abundace, 
and the nitrogen is said to be of secondary origin (Vila-Costas & 
Edmunds 1993). 

Several autors (Edmunds & Pagel 1978; Barbuy 1983; 
Tomkin & Lambert 1984; Matteucci 1986; Carbon et al. 1987; 
Henry et al. 2000) demostrated that the ratio of nitrogen to 
oxygen remains constant at lower metallicities, Z < 0.5 Z Q 
[log(0/H) + 12 < 8.3, adopting 12+log(O/H) =8.66 from 
Asplund et al. (2005)], with a plateau at log(N/0) ~ -1.5 in 
the early evolution of the galaxy, thus implying a primary ori- 
gin of nitrogen. When the oxygen abundance is greater than Z ~ 
0.5 Zo, the N/O ratio rises steeply with increasing O/H. This is 
the regime where nitrogen is predominantly secondary (Alloin 
et al. 1979; Considere et al. 2000; Pettini et al. 2008). The fact 
that the N/O ratio is relatively flat at low metallicities, indicates 
that production of nitrogen is dominated by primary processes at 
low metallicities, and by secondary processes at high metallici- 
ties (Garnett et al. 1997; Ferguson et al. 1998; Henry & Worthey 
1999). 

In order to determine the origin of nitrogen in our 1-samples, 
we estimated the nitrogen abundances for these galaxies. To 
estimate the electron temperature in the [Nn] emission region 
(^[Nn]) from log R23, we used the formula given by Thurston et 
al. (1996), 

T [m = 6065 + 1600(logfl 2 3) + 1878(logtf 23 ) + 2803(log fl 2 3), 

with r[N n ] in units of K. The ionic abundance ratio is estimated 
from the r [Nn ] temperature and the emission-line ratio ([Nn] 



AA6548, 6583)/([On] ^3727) by assuming N/O = N + /0 + , and 
the flux of [N 11] ,16548 = (1/3) [N 11] ,16583. With these, we apply 
the formula given by Pagel et al. (1992), based upon a five-level 
atom calculation: 

= to < J g" ) + 0.307 - 0.02 log f[Nll] - ™ 

where t [m = 10" 4 r [N[l] . 

In Fig. 4, the abundances of N and O are shown as a function 
of 12+log(0/H). Because nitrogen is predominantly a secondary 
element above metallicities about half solar, the galaxies of our 
samples have mainly nitrogen of secondary origin, as can be seen 
in Fig. 4. This is because we are working with both massive, 
log(M JMr /M G ) > 10.5 (see Fig. 6), and high metallicity galaxies. 
In spite of our high metallicities, we can observe signs of a hor- 
izontal population in 1-sample a at log(N/0) ~ -1.2, in a range 
of metallicities from 8.4 to ~8.6, which is a little higher than 
the standard value [log(N/0) ~ -1.5], but comparable to the re- 
sults obtained by Pettini et al. (2008) and Considere et al. (2002), 
who shown a population of galaxies in the same region. This re- 
sult is not surprising, because, as argued above, the transition 
zone from primary to secondary production of nitrogen occurs 
at 12+log(0/H) ~ 8.3, which is close to our lower limit. 

Because low metallicity galaxies (12+log(0/H) < 8.3) are 
absent in our samples, nothing can be concluded for the primary 
production of nitrogen. 

5. Discussion 

To investigate whether the origin of the decrement in metallicity 
is due to instrumental effects, to an inherent property of the sam- 
ple of galaxies, or a mixture of both, it is necessary to explore 
two important effects: the luminosity and mass-metallicity rela- 
tions, and the effect of the 3 arcsec diameter of the Sloan fibers. 

5.1. Effect of the Mass and Luminosity-Metallicity relations 

The luminosity-metallcity (L - Z) relation was first observed by 
McClure & van den Bergh (1968) in elliptical galaxies, and con- 
firmed by Garnett & Shields (1987), while the mass-metallicity 
relation was first identified for irregular and blue compact galax- 
ies by Lequeux et al. (1979), and Kinman & Davidson (1981), 
respectively, and confirmed by Skillman et al. (1989). Since 
then, as luminosities are easier to estimate than masses, many 
studies have focussed on the L - Z relation (e.g., Skillman et al. 
1989; Brodie & Huchra 1991; Zaritsky et al. 1994; Garnett et 
al. 1997; Lamareille et al. 2004, 2006; Maier et al. 2004), which 
correlates the absolute magnitude of galaxies with metallicity, 
being the more metal rich, the more luminous. The M - Z and 
L-Z relations have been studied at both, low and high redshift 
(e.g. T04; Savaglio et al. 2005; Erb et al. 2006; Maiolino et al. 
2008; Lamareille et al. 2009). 

There are two main ways to explain the origin of the M - Z 
relation, one is attributed to metal and baryon loss due to gas 
outflow, where low-mass galaxies eject large amounts of metal- 
enriched gas by supernovae winds before high metallicities are 
reached, while massive galaxies have deeper gravitational poten- 
tials which helps retaining their gas, thus reaching higher metal- 
licities (Larson 1974; Dekel & Silk 1986; MacLow & Ferrara 
1999; Maier et. al. 2004; T04; De Lucia et al. 2004; Kobayashi 
et al. 2007; Finlantor & Dave 2008). A second scenario to ex- 
plain the M - Z relation is by assuming low star formation ef- 
ficiencies in low-mass galaxies (Efstathiou 2000; Brooks et al. 
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8.97 


9.03 
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8.92 


8.92 
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Table 3. Mean and mode of the metallicity distributions, and the metallicity range that include the 66% of the total metallicity 
distribution around the mode for the different luminosity samples. 




Fig. 4. Abundances of N and O ratios as function of their 12+log(0/H) abundances derived from the Rn calibration for our three 
1-samples in redshift, following the same code of symbols used in Fig. 2. The dot-dashed line and the dashed line represent an 
aproximation of empirical limits of the primary and secondary levels, respectively, of N production, taken from Vila-Costas & 
Edmunds (1993) 



2007; Mouhcine et al. 2008; Tassis et al. 2008; Scannapieco et 
al. 2008). 

As pointed out in the high-resolution simulations of Brooks 
et al. (2007), supernovae feedback plays a crucial role in lower- 
ing the star formation efficiency in low-mass galaxies. Without 
energy injection from supernovae to regulate the star formation, 
gas that remains in galaxies rapidly cools, forms stars, and in- 
creases their metallicity too early, producing a M - Z relation 
too flat compared to observations. However, Calura et al. (2009) 
reproduced the M - Z relation with chemical evolution models 
for ellipticals, spirals and irregular galaxies, by means of an in- 
creasing efficiency of star formation with mass in galaxies of all 
morphological types, without the need of outflows favoring the 
loss of metals in the less massive galaxies. A recent study that 
supports this result for massive galaxies, is the one of Vale Asari 
et al. (2009), modelling the time evolution of stellar metallicity 
using a closed-box chemical evolution model. They suggest that 
the M - Z relation for galaxies in the mass range from 10 9 8 to 
j Qii.65 jy[ Q j s mam iy driven by the star formation history and not 
by inflows or outflows. 

As explained in Sec. 2, we selected small intervals of lu- 
minosity for all redshift-samples. Such selection was aimed to 
avoid possible biases since the SLOAN star-forming sample be- 
comes incomplete at redshifts above z > 0.1 (e.g., Kewley et al. 
2006). Then, this could introduce a bias, since our high-redshift 
sample is formed by the most luminous galaxies, thus resulting 
in higher metallicity estimates. 

As can be seen from Fig. 5, the galaxies of our samples do 
not show a luminosity-metallicity dependence, since this rela- 
tion can be clearly seen only when spanning ranges of more 



than ~4 magnitudes in luminosity (e.g. T04). The masses of 
our galaxies were estimated in order to check their behaviour 
in our luminosity intervals with the STARLIGHT code, using a 
Chabrier (2003) initial mass function between 0.1 and 100 M Q , 
for details on the mass estimates, see Mateus et al (2006). The 
M - Z diagram for our 1-samples do not show any correlation 
either as can be seen in Fig. 6, since masses again correspond to 
a small interval in luminosity for each redshift bin. 

For each 1-sample, we overplot the polynomial fit of the lo- 
cal M - Z relation of Tremonti et al. (2004), which is valid over 
the range 8.5 < log(M.„ f „/M ) < 11.5, with a steep M - Z re- 
lation for masses from 10 8 ' 5 to 10 10 5 M , that flattens at higher 
masses. As can be seen in Fig. 6, galaxies of the 1-sample a, b 
and c, correspond to this flat and massive zone [log(M JtoJ /M ) 
> 10.5] of the M - Z relation, for the reasons discussed above. 
The mass ranges were we can found two thirds (~66%) of each 
1-sample are shown in Table 4, as well as the median Petrosian r 
magnitude, and the median mass in log(M s/flf /M ). 

5.2. Effects of the Sloan fiber diameter 

With respect to the Sloan fiber diameter, and depending on the 
galaxy size, we expect that at higher redshift the projected size 
of the Sloan fibers will cover a larger fraction of galaxy area than 
for nearby galaxies. This effect, as argued in Sect. 3, could in- 
troduce a bias in our samples since integrated metallicities are 
lower than nuclear ones. In order to quantify this contribution, 
we estimate the percentage of angular size of each galaxy inside 
the three arcsec diameter of the Sloan fiber using the Petrosian 
total radius in r-band in arcsec. To this aim, we divide the fi- 
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Table 4. Median Petrosian r magnitude (M r ), median of the logarithm-mass [log(M SM ,7M )], and the logarithm-mass range that 
include the 66% of the total mass distribution around the median logarithm-mass value for 1-samples a, b, and c. 




Fig. 5. Absolute Petrosian r magnitude and 12+log(0/H) for our 1-sample of galaxies (L - Z relation). Symbols follow the same 
code used in Fig. 2. 




Fig. 6. Relation between stellar mass, in units of solar masses, and 12+log(0/H) for our 1-sample of galaxies (M - Z relation). 
Dashed line represent the fit of Tremonti et al. (2004), and symbols follow the same code used in Fig. 2. 



bre radius (1.5 arcsec) by the Petrosian total radius in r-band 
(petrorad r ), as can be seen in the histograms of Fig. 7 for all 
our 1-samples. Thus, this ratio can be taken as the fraction of the 
galaxy size that is actually covered by the fibre. 

As expected, in all 1-samples the distribution shifts to a maxi- 
mum coverage of the galaxy size as redshift increases. This is an 
effect that must be taken into account, because it could change 
the metallicity estimation according to the fraction of galaxy di- 
ameter inside the Sloan fiber. As argued in L09, the decrement 
observed in metallicity for 1-sample c can not be atributed to the 



percentage of galaxy inside the Sloan fiber, because fibers cover 
less than 50% of galaxy sizes for ~95% of this sample. However, 
the effect could be noticeable at lower redshifts. 

To minimize systematic errors from this aperture biasses, as 
explained in Sec. 2, we selected galaxies with z > 0.04 follow- 
ing the recommendations of Kewley et al. (2005), who inves- 
tigate the effects of a fixed-size aperture on spectral properties 
for a large range of galaxy types and luminosities, concluding 
that a minimum flux covering fraction of 20%, corresponding 
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Fig. 7. Percentage of angular size of each galaxy inside the three arcsec diameter of the Sloan fiber for 1-samples a, b, and c. 



to a median redshift of z ~ 0.04, is required for metallicities to 
approximate the global values. 

It is important to point out that in Fig. 7 we are plotting the 
percentage of angular size of each galaxy inside the three arcsec 
diameter of the Sloan fiber, which is not the same as the flux 
covering fraction used by Kewley et al. (2005). 

However, in spite of this redshift limit and in order to test 
how much the percentage of angular size inside the Sloan fiber 
diameter affect our metallicities, we compare in Fig. 8 our orig- 
inal metallicity estimates of redshift-sample zo, with that from 
galaxies with a percentage of angular size within the fiber >20%, 
which corresponds to a 45% of the sample. We observe that the 
result is the same, with a quite small difference in the mean of the 
order of ~0.001 in 12+log(0/H) for galaxies with a percentage 
of angular size within the fiber >20%. 

If redshift sample zo, which contains the maximum fraction 
of galaxies with a percentage of angular size <20% inside the 
Sloan fiber diameter, do not show variations in the mean metal- 
licity, we can assume that none of our samples are significantly 
affected by aperture effects due to the Sloan fiber diameter. 
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Fig. 8. Normalized metallicity histograms for redshift-sample zo- 
Dotted clear line is the same as in Fig. 3a, and dotted dark line 
are galaxies of the same redshift-sample with covering fraction 
>20. 



6. Conclusions 

Although similar studies of metallicity using SDSS exists, they 
are either restricted to a redshift ~0. 1 (e.g. T04), or do not seg- 
regate their samples as a function of redshift (e.g. L06), thus 
making not possible to detect metallicity evolution at low red- 
shift. On the other hand, studies at high redshift are statistically 
limited. 

We divided our sample in redshift intervals of Az ~0.1 with 
the goal of identifying any evolution in metallicity, and each red- 
shift in small intervals of luminosity in order to avoid biases due 
to the luminosity and mass-metallicity relations. Because we are 
comparing galaxies from redshift 0.04 to 0.4, we selected lumi- 



nous and massive galaxies [~ log(M.„ flr /M ) > 10.5], which are 
present in all redshift intervals. 

We conclude that the nitrogen production for our sample of 
galaxies is mainly secondary, because our sample is formed by 
massive, luminous and high metallicity galaxies. We remark that 
in this work we can not conclude nothing about the production 
of nitrogen in low massive galaxies. 

We derived the M - Z and L - Z relations for our sample of 
galaxies, showing a flat L-Z relation, and a M - Z relation pop- 
ulated only in the massive zone ~ log(M sfflr /M ) > 10.5, which 
is due to the absence of the dwarf galaxies. Both relations be- 
have as expected since our sample is selected to cover the more 
luminous and hence massive galaxies. 
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Since we have metallicity estimates for redshifts up to 0.4 in 
bins of 0.1, we are able to represent redshift versus metallicity, 
as shown in Fig. 9. We plotted the mean metallicity for all 1- 
samples a, b, and c. Error estimates of the mean metallicity were 
obtained from the line error fluxes provided by STARLIGHT, 
taking into account the number of galaxies of each sample. 

As shown in Fig. 9, we can observe, for 1-sample a, an ini- 
tial metallicity of ~9.02 for zo, followed by a small decrease in 
metallicity for z\- For 1-sample b, metallicity of z\ and zi re- 
mains constant, a trend also observed for 1-sample c in the same 
redshift-samples, but with lower metallicity values, and finally 
an important decrement for Z3- We overplot the model of Buat 
et al. (2008) for galaxies with a rotational velocity of 360 km/s, 
which correspond to a log(M Jfa ,/M ) ~ 1 1.2 in his model. 

Since the model of Buat et al. (2008) is calibrated in solar 
metallicities, with log(Z/Z ) ~ -0.03 at redshift zero, we assume 
a solar metallicity of 12+log(O/H)=9.05 in order to match his fit 
to our metallicity values, finding a good correlation between her 
model and our metallicity dependence versus redshift. Finally, 
we fit a second order polynomial (see Fig. 9). 

It is important to have in mind that this result is valid for 
massive, luminous, and high metallicity galaxies. Lower metal- 
licity and stellar mass galaxies are absent from our sample due 
to the selection criteria applied. 

Although it is well known that metallicities decrease with 
redshift, it is the first time that a statistically significant sample 
of galaxies is analized looking for metallicity evolution at such 
a low redshift, detecting small decrements as redshift increase, 
with prima facie evidence of a quite important decrement at the 
redshift-interval 0.3-0.4, as already pointed out in L09. 



9.05 r 




8.90 - 



0.0 0.1 0.2 0.3 0.4 

redshift 

Fig. 9. Metallicity evolution derived from our 1-samples up to 
redshift 0.4, vertical lines show the error with respect to the mean 
metallicity for every redshift interval. Circles, triangles and as- 
terisks, represent the mean metallicity of the 1-samples a, b and 
c, with its respective redshift-samples. Solid line represents the 
model of Buat et al. 2008 for galaxies with a rotational velocity 
of 360 km/s, and dashed line represents a second-order polyno- 
mial fit (y=ao + a\x + a2X 2 ), with ao = 9.018, a\ = 0.015, and 
a 3 = -0.799. 
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